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Abstract 


To develop advanced control systems for optimiz- 
ing aircraft engine performance, unmeasurable output 
variables must be estimated. The estimation has to be 
done in an uncertain environment and be adaptable to 
varying degrees of modeling errors and other varia- 
tions in engine behavior over its operational life cycle. 
This paper presents an approach to estimate unmea- 
sured output variables by explicitly modeling the ef- 
fects of off-nominal engine behavior as biases on the 
measurable output variables. A state variable model 
accommodating off-nominal behavior is developed for 
the engine, and Kalman filter concepts are used to es- 
timate the required variables. Results arc presented 
from nonlinear engine simulation studies as well as the 
application of the estimation algorithm on actual flight 
data. The formulation presented has a wide range of 
application since it is not restricted or tailored to the 
particular application described in the paper. 

Nomenclature 


A,B,C, Dy 

F, Gy H 

Aj 

b 

CDF 


system matrices in state 
variable representation 

nozzle area, in. 
bias vector 

component deviation factors 
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CIVV 

DEEC 

DmOZ 

Dram 

d 

E 

EMD 

e 

Fa 

Fmp 

I 

K 

Ni 

N 2 

P 

PD 

PLA 

PSC 

PSM 

Pt 2 5 

Pt s 

Pn 

PTi 

Q 

R 

Rcvv 

SM f 


fan inlet guide vane angle, deg 
digital electronic engine control 
nozzle drag, lb 
ram drag, lb 
difference 

expectation operator 
engine model derivative 
state error vector 
gross thrust, lb 

net propulsive thrust, lb 
identity matrix 

Kalman filter gain 
fan rotor speed, rpm 
core rotor speed, rpm 
Riccati matrix 
burner static pressure, lb/in 2 
power lever angle, deg 
performance seeking control 
propulsion system model 
compressor inlet total pressure, lb/in 2 
low turbine inlet pressure, lb/in 2 
afterburner inlet total pressure, lb/in 2 
nozzle throat total pressure, lb/in 2 
state noise covariance matrix 
measurement noise covariance matrix 
compressor stator vane angle, deg 
fan stall margin 



SMhc 

high compressor stall margin 

SOAPP 

state-of-the-art propulsion program 

SSM 

steady state model 

SVM 

state variable model 

TMT 

turbine metal temperature, °R 

^2 5 

compressor inlet total temperature, °R 

Tt> 

burner inlet total temperature, °R 

Tn 

burner exit total temperature, °R 

T T<5 

low turbine inlet total temperature, °R 

Tt 6 

afterburner inlet total temperature, °R 

t Ti 

nozzle throat total temperature, °R 

u 

control input vector 

WCfan 

corrected fan air flow, lb/sec 

WChpc 

corrected compressor air flow, lb/sec 

W F 

main burner fuel flow, lb/hr 

W\ 

state excitation noise 

W2 

measurement noise 

X 

state vector 

y 

output vector 

Vaux 

vector of auxiliary, (unmeasured) 
output variables 

Z 

augmented state vector 

8 

variation from trim values 

Superscripts 

T.' 

transpose of a matrix 


estimated value of variable 


derivative 

Subscript 

b 

output bias term 

m 

flight or simulated data 

N.L. 

nonlinear 

i 

trim, initial, or steady state 

1 

augmented system matrices 


Introduction 

Efforts to improve aircraft turbine engine efficiency 
have led to an increase in the number of engine control 


variables and a corresponding increase in the complex- 
ity of control laws. Control laws for current engines 
are based on classical control theory and empirical 
schedules for a nominal engine. Classical control the- 
ory has served well for the current and older engines. 

The design of future fighters as multifunction aircraft 
and development of integrated flight/propulsion con- 
trol systems, however, require sophisticated control 
systems capable of obtaining the maximum perform- = 
ance from the engine. Optimal control techniques us- 
ing modem control theory are required to obtain ad* 
ditional gains in engine performance. For modem air- 
craft, accounting for engine variations through designs * 
based on predetermined control schedules is increas- 
ingly difficult because of the increased complexity and { 
increased number of control effectors on the engines. f 
Engine-to-engine component variations, engine dete- i 
rioration, and off-nominal behavior are difficult to ac- 
count for in the design of control system schedules. 

An adaptive control algorithm, which computes op- = 
timal control trim settings for the engine while maxi- j 
mizing the vehicle performance for a given flight con- 1 
dition, accounts for these variations better than gain 
scheduling. Specifically, an adaptive trim control sys- 
tem computes and applies an incremental steady state 
trim to enhance the engine performance. 1 

For over a decade, the National Aeronautics and 
Space Administration (NASA) Ames Research Cen- 
ter, Drydcn Flight Research Facility (Ames-Dryden) 
has conducted a multidisciplinary flight research pro- 
gram on an F- 15 airplane. Significant portions of this 
research involved the flight evaluation of advanced 
propulsion control concepts in programs such as digi- 
tal electronic engine control (DEEC), the FI 00 engine 
model derivative (EMD), and highly integrated digi- 
tal electronic control (HIDEC). 2 The increased perfor- 
mance and improved fuel economy demonstrated on 
the F-15 HIDEC research vehicle is the basis of the 
performance seeking control (PSC) program, which 
will provide additional improvements in these areas. 

Ames-Dryden, McDonnell Aircraft Company, and 
Pratt & Whitney are currently developing and demon- 
strating an adaptive PSC system in flight on a NASA 
F-15 airplane powered by FI 00 EMD engines. The 
PSC system optimizes aircraft performance by apply- 
ing adaptive trim control to the propulsion system op- - 
crating in a pseudo-steady-state cruise mode. The 
trim schedules are determined for a highly nonlin- 
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ear propulsion system which has system and measure- 
ment noise, unmeasurable parameters, and sensitivity 
to normal deterioration over its life cycle. 

Figure 1 shows the adaptive trim control structure 
used for the PSC. The state variable model (SVM) 
and the steady-state model (SSM) which model the dy- 
namic and steady-state behavior of a nominal engine, 
arc key components of the system. These models are 
stored onboard the aircraft in a table look-up form and 
are discussed in more detail in the following section. 

These models are used in formulating the propulsion 
system model (PSM) which represents a small pertur- 
bation model of the actual (light propulsion system. 
The PSM contains relations which provide estimates 
of performance measures (such as augmentor effects, 
thrust, and stall margins) and constraint equations. A 
linear programming algorithm is used to find the opti- 
mal solution and these commands are then applied to 
the engine through the DEEC. 

The values of output variables, which are often not 
directly measurable, are needed for the optimization 
algorithm used in the PSC. These variables are esti- 
mated under changing levels of engine health, man- 
ufacturing differences between engines, and other off- 
nominal behavior. Accommodating these performance 
variations in engines has been investigated in two re- 
cent studies. 3,4 

Reference 3 presents an algorithm for estimating the 
cause and level of off-nominal engine operation by us- 
ing a Kalman filter algorithm to estimate five engine 
factors. These five factors, referred to as component 
deviation factors (CDF), compensate for off-nominal 
performance. These factors were estimated by treat- 
ing them as biases, and the original state vector was 
augmented to give five additional states. 5 These five 
factors are not explicitly used in the optimization algo- 
rithm and their physical significance is unclear because 
the formulation does not account for biases, prediction 
errors, and Reynolds number effects. Since the coeffi- 
cients with respect to the CDF parameters arc required 
in the Kalman filter development, the CDF formulation 
requires detailed modeling of the off-nominal process. 
A flight data evaluation of this algorithm is described 
in Ref. 6. 

In Ref. 4, a component tracking filter is used to 
achieve the model accuracy required to optimize en- 
gine performance. The component tracking filter com- 


bines the concept of state tracking and adaptive filter- 
ing to minimize engine/model mismatch. It is based on 
a frequency decomposition of the differences between 
the sensed engine parameters and the model values. 

This paper presents another method of accounting 
for off-nominal operation and other modeling inaccu- 
racies. Since any variation from the nominal model 
would result in a change in the sensed values of the 
measured outputs, the off-nominal behavior of an en- 
gine is characterized in terms of these changes. Uncer- 
tainties associated with any given engine will be repre- 
sented as systematic errors in the sensed output param- 
eters. These systematic errors will be accounted for by 
augmenting the original state equation with bias states. 
A Kalman filter is used to estimate the original engine 
states and the bias states. The Kalman filter inputs are 
measurements from standard FI 00 engine control in- 
strumentation. The auxiliary output equations for the 
unmeasured output variables are modified to include 
the effect of the bias states. 

The concept is validated by applying the developed 
filter on both simulation and flight data. For the sim- 
ulation data case, the output variables were estimated 
by using the data from the available nonlinear engine 
simulation. Both a nominal engine and an engine in 
which intentional degradation was introduced to create 
off-nominal behavior were considered. For the flight 
data case, the estimation process was performed using 
actual flight data from an F-15 aircraft. For this case, 
comparative results are also presented for the proposed 
algorithm and the CDF formulation. Both the sim- 
ulation and flight evaluations were carried out for a 
flight condition of Mach 0.90 and 30,000 ft, for a part 
power setting. 

Engine Description 

The engine used in this study is the Pratt & Whitney 
F100 EMD low-bypass ratio, twin spool, afterburn- 
ing turbofan engine 7 (Fig. 2). The engine is controlled 
by a DEEC, a full-authority digital electronic control 
system which performs the functions of the standard 
F100 engine hydromechanical, unified fuel control, 
and supervisory digital electronic engine control. 

Engine Models 

Pratt & Whitney has developed a comprehensive 
nonlinear dynamic engine model, the state-of-the-art 
propulsion program (SOAPP) model. This model is 


the best representation of the engine and predicts en- 
gine performance with minimal error over the full 
power range and (light envelope and for both steady- 
state and transient operation. This nonlinear simula- 
tion is a high-fidelity model that represents each com- 
ponent in the engine and control but does not run in 
real time. 

For real-time use, a set of linearized SVMs were de- 
veloped from the SOAPP model. To cover the entire 
flight envelope, 49 models were developed. The model 
is selected as a function of burner static pressure (PB). 
These models compare well with the large scale non- 
linear aerothermal model and actual engine test data, 
and they can be implemented efficiently in real time. 
Figure 3 shows a simulation model for the FI 00 engine 
based on the state variable formulation. 

The SSM engine relationships and trim predictions 
(basepoints) are also derived from the SOAPP model. 
A two-dimensional table look-up scheduled on 7 val- 
ues of PB and 40 values of afterburner total pressure 
(Pt 6 ) is needed to represent the steady state informa- 
tion. Each SSM consists of a basepoint control vector, 
a basepoint output vector, and a sensitivity coefficient 
matrix which relates the changes in control positions 
to change in outputs. 

The PSC algorithm requires the variables listed in 
Table 1, which are functions of the engine states and 
the input control variables. These variables include 
engine outputs which cannot be measured but are re- 
quired to calculate performance measures of the en- 
gine. An additional set of variables, which are non- 
linear functions of the unmeasured output variables, 
arc listed in Table 2. These variables are used to pre- 
dict both the engine performance and the constraints 
needed to develop optimal engine controllers. 

Kalman Filter Concepts 

The entire state vector of the system to be controlled 
is often assumed to be measurable. Most of the so- 
lutions to optimal control problems are obtained as a 
feedback law implcmentable only if the entire state 
vector is available. In most complex systems the en- 
tire state vector cannot be measured, and a suitable ap- 
proximation to the state vector must be determined and 
substituted into the control law. The system that pro- 
duces, in deterministic setting, an approximation to the 
state vector is called an observer. 8 


Kalman and Bucy solved the optimal observer prob- 
lem in a stochastic environment, and this solution has 
had a tremendous impact on optimal filtering theory. 9 
The Kalman filter represents the most widely applied 
and demonstrably useful result to emerge from the 
state variable approach of “modem control theory.” 10 

The system is 

x - Ax + Bu -f vu\ (1) 

y = Cx + Du + w 2 (2) 

Where A , B , C, and D are system matrices in state 
variable representation, x is the state vector, u is the 
control input vector, y is the output vector, w\ is the 
state excitation noise, and W 2 is the observation or 
measurement noise. Both w\ and wi are white, un- 
correlatcd Gaussian processes, with intensity Q and 
R respectively. 

The observer is 

x - Ax + Bu + K[y — Cx — Du] 

where K is the Kalman filter gain. 

The optimal observer problem is finding the matrix 
K so as to minimize E{e J Re}, where 

e = x — x 

and R is a positive-definite symmetric weighting ma- 
trix. In this problem, E is the expectation operator and 
e is the state error vector. If R is a positive-definite 
matrix, the optimal observer is called nonsingular. The 
Kalman filter is the solution to the nonsingular optimal 
observer previously outlined. The optimal observer 
problem is solved by choosing the gain matrix. 11 

K = PC 7 R~ l 

where P is the state error covariance matrix, 
E[(x - x)( x — x) J ], and is the solution to the matrix 
Riccati equation 

P = AP + PA J + Q - PC T R~ l CP 

For a time invariant case, the steady state solution for 
Pisa constant matrix and is a unique nonnegative def- 
inite solution of the algebraic Riccati equation 

0= AP+ PA J + Q - PC T R~ l CP 
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Figure 4 shows a typical Kalman filter structure used 
to estimate states and outputs. 

Proposed Formulation 


where P is the steady state solution to the 
Riccati equation 

0 =AiP+ PA] + GQG J - PCjir'CiP 


In Kalman filter derivation, linear models for the 
system dynamics and measurement relation are as- 
sumed to be adequate for developing optimal estima- 
tors. No model is perfect, and a linear model, in 
particular, is the result either of intentional approxi- 
mation and simplification or of a lack of knowledge 
about the system being modeled. 12 To account for de- 
graded engine operation and modeling inaccuracies, 
the proposed formulation augments the output vec- 
tor by adding a bias vector to represent the uncertain 
parameters. 5 The dynamic equations can thus be ex- 
pressed as 


x = Ax + Bu + w\ 
y = Cx + Du + 6 + W2 


where b is the bias vector. The bias vector is estimated 
by adjoining 6 to x and defining a new state vector, 2 


x 


2 = 


b 


with the condition 

6 = 0 


The state equation can be rewritten as 

2 = A\z + B\u + Gw 1 
y = C\z + Du + W 2 


where 


Ai = 


A 

0 


0 

0 


B\ = 


G = 


B 

0 

I 


Cl = [C : I ] 

If the estimate of 2 is 2 , where 

x 
i 

2 = 


then the Kalman filter estimate is given by 

z = A\z + Diu+ PC] R-'lv-Ciz - Du] 


The auxiliary set of unmeasured output variables 
(y ou *) arc related to the engine states and control in- 
puts through the algebraic equation 

Vaux = Hz + Fu 

Details of the state variable formulation for the 
F100 engine are presented in the appendix. The (y 0UI ) 
outputs arc listed in Table 1 . 

In spite of the mathematical formalism of the 
Kalman filler, engineering insight and experience is 
required to develop an effective operational filter al- 
gorithm. A mathematical model of both the system 
structure and uncertainty is inherently embodied in the 
Kalman filter structure. The main design problem is at- 
taining an adequate mathematical model upon which 
to base the filler. Even after selecting an appropriate 
model, the matrices Q and R can be difficult to de- 
termine. This is done by a process called “tuning” the 
Kalman filter. It is a trial and error procedure for deter- 
mining which matrix values yield the best estimation 
performance for that particular filter structure. 

The matrix R was determined by analysis of (light 
data available for the F100 engine. The elements of 
matrix Q were, however, selected by evaluating the 
performance of the Kalman filter by trial and error. 
Figure 5 shows the implementation process used to es- 
timate the output variables for the FI 00 engine using 
the Kalman filter. 

This proposed formulation estimates unmeasured 
output variables by explicitly modeling the effects of 
off-nominal engine behavior as biases on the measur- 
able output variables. 


Results 

The proposed estimation algorithm was developed 
and evaluated for a Mach 0.90 and 30,000 ft flight con- 
dition. The algorithm was evaluated by a comparison 
with SOAPP simulation results and also by application 
to flight data. The flight data results were compared 
with the CDF formulation results for the same data. 

Simulation Evaluation 

The SOAPP simulation evaluations consisted of es- 
timating the desired variables using both a nominal and 
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a degraded engine. In each case, the power lever angle 
(PL A) was held to 37° for 15 sec and then stepped up 
to 43° and held constant for the remainder of the run. 

Measured outputs were obtained from the SOAPP 
simulation and were corrupted with noise, as shown 
in Table 3. These are typical values obtained from 
flight data. The measurements with noise and the 
values of the control variables were entered into the 
estimation algorithm and the desired estimates were 
obtained. The Kalman filter state vector, a perturba- 
tion of the steady state conditions, was initialized to 
zero for all states. 

The algorithm needed to generate consistent state 
estimates which were robust with respect to the mea- 
surement covariance matrix Q (the only variable se- 
lected by trial and error). An important aspect of the 
development is determining unmeasured output vec- 
tor, £oux* Inconsistent estimates of the states would 
give different values of y aux for different values of Q 
when applied to the same data. 

The state vector estimates converged to the same 
value for different values of Q. This was evaluated 
for values of Q = I and Q = 101. The difference in the 
estimated states for Q = I and Q = 101, for a nom- 
inal engine, is shown in Fig. 6. This figure shows 
that the state estimates converge to the same value and 
the effect of change in Q on the steady-state response 
is minimal. 

The five measured output variables obtained from 
the SOAPP for a nominal engine were compared with 
the estimates of these variables obtained from the filter 
(Fig. 7(a)). The prediction values subtracted from the 
simulated measurements were held constant through- 
out the run. These values were the same as the sim- 
ulated measurements at the beginning of the run, ac- 
counting for the excellent comparison over the initial 
interval. The Kalman filter was not updated in this 
evaluation, so the comparisons indicate that the model 
is quite robust. The comparisons are very good in spite 
of the large change in the operating conditions. The 
CDF based formulation would have used five different 
models for the PD change of this maneuver. 

Figure 7(b) shows the measurement bias estimates. 
As expected, they are nearly zero until the PLA is in- 
creased. As the engine attains a new operating condi- 
tion, the bias parameters increase to levels which ac- 
count for the effects not modeled in the SVM. 


To assess the condition when significant differences 
exist between the measured data and the predicted 
data, the following nominal biases were added to 
the simulated flight data: A Ah (fan rotor speed) = 
50.0, A Ah (core rotor speed) = 50.0, APS = 2.0, 
A Tt 4 5 (low turbine inlet total temperature) = 30.0, and 
A Pp 6 = 0.5. The results of this evaluation (Fig. 8(a)) 
show that the tracking of the five measurements is 
again very good. The final values of the bias estimates 
(Fig. 8(b)) are the sum of biases estimated in Fig. 7(b) 
and the biases placed on the simulated measurements 
as previously listed. 

In Fig. 9, estimates of the unmeasured output vari- 
ables (y flUX ) are compared with the actual values ob- 
tained from the SOAPP. The estimates show good 
tracking of the simulation values. 

Simulation evaluations were then carried out for 
a degraded engine by simultaneously introducing the 
following deteriorations: (a) high turbine efficiency is 
2.5 percent below nominal, (b) low turbine efficiency 
is 2.5 percent below nominal, (c) compressor airflow 
deviation is 1 lb/sec less than nominal, and (d) the fan 
airflow deviation is 5 lb/sec less than nominal. 

The results for the simulated degraded engine arc 
presented in Fig. 10. These results are similar to the 
results of Fig. 7 and demonstrate the adaptability and 
robustness of the proposed estimator to degraded en- 
gine performance. Again, the Kalman filter was not 
updated during the evaluation and the predicted con- 
stant values subtracted from the simulated data were 
the same as those for an engine that was not degraded. 

Flight Data Evaluation 

The Kalman filter formulation was also evaluated on 
flight data obtained on the NASA F-15 research air- 
craft. The flight data was obtained at Mach 0.90, an 
altitude of 30,000 ft, and a PLA of 43.5°. The time 
history of the test data (Fig. 1 1) starts with no bleed air 
being extracted from the test engine. Approximately 
40 sec into the run, the pilot manually changed the 
bleed switch to extract all the aircraft bleed air require- 
ments from the test engine. This maneuver was de- 
signed to simulate a change in engine operating effi- 
ciency. The engine control system increased fuel flow 
{Wp) to maintain the scheduled fan speed, resulting in 
an increase in Tt a 5 ■ After holding this bleed condition 
for approximately 70 sec, the bleed was again switched 
back to the initial no bleed air condition. 
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The Kalman filter estimation results are shown in 
Fig. 12. Figure 12(a) shows that the filter tracks the 
flight measurements accurately. Initial discrepancies 
occur because the bias estimates start at zero; how- 
ever, this startup transient is brief, with good track- 
ing occurring in approximately 20 sec. Although 
the tracking quality is slightly worse at the time the 
bleed switching occurs, the filter rapidly adapts to 
the simulated change in engine efficiency. The bias 
estimates, shown in Fig. 12(b), converge rapidly to 
steady-state values as the engine state is changed from 
one condition to another. The initial startup tran- 
sient could be minimized by initializing the bias esti- 
mates with the actual values of the biases for the given 
flight condition. 

Figure 13 shows the results from the proposed 
formulation compared with the corresponding results 
from the CDF formulation. The results were obtained 
using the flight data shown in Fig. 11. The results show 
that the performance obtained by the proposed method 
compares favorably with the CDF procedure. A signif- 
icantly improved startup transient performance is evi- 
dent. Figure 14 presents similar comparisons for the 
estimates of normally unmeasured output variables. 
Figure 14(a) shows the estimate of compressor inlet to- 
tal temperature ( Tt 25 ) and the measured values. The 
superiority of the proposed formulation is clearly ev- 
ident, if the measurement of Tt 2 5 is considered reli- 
able. Figure 14(b) shows the comparative estimates of 
corrected fan airflow ( WCfan )• The values are com- 
parable, with better transient performance for the pro- 
posed formulation. 

Concluding Remarks 

An approach has been proposed to estimate the un- 
measured or auxiliary output variables of a turbofan 
F100 engine by using Kalman filter concepts. The 
approach is based on explicitly modeling the effects 
of off-nominal engine behavior as biases on the mea- 
sured output variables. Results are presented for esti- 
mates of the output variables and are compared with 
values obtained from detailed nonlinear simulation of 
the engine. The evaluation was carried out for both 
a nominal engine and an engine in which intentional 
deterioration was introduced. The proposed filter was 
also evaluated for output estimation using actual F-15 
flight data. 


The formulation is robust with respect to the value 
of state covariance matrix Q. A critical component 
of the performance seeking control (PSC) problem for 
the FI 00 engine is determining consistent values for 
auxiliary output variables. Consistent estimates for the 
states were obtained for different values of Q and thus 
consistent estimates of the auxiliary output variables 
are ensured. 

The proposed estimation algorithm was able to ac- 
curately predict the values of the output variables for 
the simulation studies for both nominal and degraded 
engine conditions. The proposed algorithm has been 
validated by comparing its estimates with the values 
from the detailed nonlinear simulation, and it has per- 
formed well on flight data. A comparative study of 
the proposed algorithm results with component devia- 
tion factors (CDF) results gave additional proof of the 
validity of the concept. Unlike the CDF method, the 
proposed algorithm does not require detailed model- 
ing of the engine degradation process. This formula- 
tion has a wide range of application because it is not 
restricted or tailored to the particular application de- 
scribed in this paper. 
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Appendix — State Variable Auxiliary Output 

For the system being considered, the complete state 
variable model is 

8z = A\8z + £?i< 5 u + Gvj\ 

8 y = C\ 8z + D8u + ui 2 


where 1 indicates augmented system matrices, and 



8N ! 





5z = 

8N2 

STMT 

N 2b 

Su = 

8W f 

8Aj 

8CIVV 

8y = 

8N 1 

6 N 2 

8Pn 

8PB 

ST n , 


Pt H 

PB 


8 RCVV 









where TMT is the turbine metal temperature, b de- 
notes the output bias term, CIVV is the fan inlet guide 
vane angle, RCVV is the compressor stator vane an- 
gle, and 


Estimation Formulation for an F100 Engine 

where A, B, C, and D are constant perturbation ma- 
trices, numerically derived from the SOAPP, w\ is the 
state noise with covariance Q, and w 2 is the measure- 
ment noise with covariance R. The elements of R are 
obtained from a priori flight data, while those of Q are 
selected by trial and error. 

The auxiliary set of unmeasured output variables 
iVaux) listed in Table 1 , is given by 

I/oui — H8z -t- F 8u "t yi 

where 

H = [Hi ; H 2 ] 

and H 2 reflects the effect of estimated biases and its 
elements arc derived from the SVM, H\ and F are 
perturbation matrices derived from the SOAPP, and y, 
is the vector of predicted trim values for the auxiliary 
output variables, which is obtained from the SVM. 


A\ 


A : 0 

0 ; 0 


a = [c ; i] 


B 1 = 


G = 


B 

0 

I 

0 
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Tabic 2. Nonlinear engine variables. 


Table 1. Linear auxiliary output variables, PSC 
algorithm requirements. 


^2,5 

compressor inlet total pressure 

PB 

burner static pressure 

Pn 

afterburner inlet total pressure 

^2 5 

compressor inlet total temperature 

Tt 3 

burner inlet total temperature 

Tn 

burner exit total temperature 

Tn, 

low turbine inlet total temperature 

Tt 6 

afterburner inlet total temperature 

WCfan 

corrected fan air (low 

WChpc 

corrected compressor air flow 


Dmoz 

nozzle drag 

Dram 

ram drag 

Fg 

gross thrust 

Fnp 

net propulsive force 

Pr; , 

nozzle throat total pressure 

SM f 

fan stall margin 

SMrc 

high compressor stall margin 

Tt, 

nozzle throat total temperature 

Table 3. 

Measurement noise statistics. 

Parameter Standard deviation 

N X 

7 rpm 

Nz 

7 rpm 

Pn 

0.3 lb/in 

Pt* 

0.6 lbrin 

Tt< 5 

4 °R 



Fig. 1 The performance seeking control adaptive control system. 
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Fig. 4 The Kalman filter structure. 

PB ' P Tr- yv 



/ SVM trim \ 

^o b ' Pt v 

f SVM trim | 


"l tables / 


1 tables J 


Proposed 

Kalman 

filter 


Flight data | Nonlinear 
calculations. 


Z = [n,, N 2 , TMT, N, b , N 2b , P T| . b , PB b , T^J 
u= Wp, A j , C IVV , RCWj' 
y = [M,,N 2 ,PB,P T6 ,T T45 ]' 

Vaux= [ P T 2 . 5 ’ T T 2 5 . t t 3 . t t 4 . t t 6 . w Cfan. w c H pcz]’ 

N.L." [ SM F’ SM HC' D RAM' D NOZ’ F NP’ F G ,T T 7 ' P T 7 ] 

Fig. 5 Modified estimation process using the proposed Kalman filter. 
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— Measured 


-- Estimated 



N 2 > 

rpm 


11,200 

11,000 

10,800 

10,600 

10,400 







(a) Measured and estimated engine outputs. 

Fig. 8 The F100 engine simulation parameter estimates with biased measurements for a nominal engine, with 
PLA increased from 37° to 43° at 15 see. 
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— SOAPP 


— Estimated 





Fig. 9 The F100 engine simulation auxiliary output estimates for a nominal engine, with PLA increased from 37° 
to 43° at 15 see. 
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— Measured 


— Estimated 







(a) Engine output estimates. 

Fig. 10 The F100 engine simulation parameter estimates for a deteriorated engine, with PLA increased from 37° 
to 43° at 15 see. 
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Time, sec 
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(a) Measured output variables. 

Fig. 11 The F-15 airplane measured engine parameters during compressor bleed variations at Mach 0.90, an 
altitude of 30,000 ft, and PLA = 43°. 






Time, sec 


900193 


(a) Output estimates. 

Fig. 12 The FI 00 engine parameter estimates from the flight data in Fig. 11. 









(a) Compressor inlet total temperature estimates. 

Fig. 14 The proposed formulation and the CDF formulation engine parameter estimates from flight data compared 
with measured engine parameters. 
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